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Steady-state properties of 180° Bloch domain wall (DW) in superconducting ferromagnet (SCFM) 
are studied. The distribution of magnetic field above and below the surface of the SCFM due to the 
permanent magnetization supercurrent flowing in the DW plane is calculated by solving Maxwell 
equations supplemented by London equation. It is shown that part of the magnetic flux of the 
two neighboring domains closures in the nearest vicinity of the surface of the sample giving rise 
to declination of the line of the force from being parallel to the DW plane. As a result, the value 
of the normal component of magnetic field at the surface of the sample reaches only half of the 
value of the bulk magnetic flux. At the distances of the order of value of London penetration depth 
the magnetic field decreases as inverse power law due to the long-range character of dipole-dipole 
interaction. The last two circumstances are important for comparison the calculated magnetic field 
with the data obtained by the methods on measurement of normal component of magnetic field, e. g. 
Hall probe technique, aimed to confirm the existence of magnetic order parameter in unconventional 
superconductor. 

PACS numbers: 74.25.Ha, 74.90.+n, 75.60.-d 



I. INTRODUCTION 

Unconventional superconductors are a subject of in- 
tensive experimental and theoretical studies during past 
decade [lf-@- Coexistence of superconductivity and 
magnetism results in a number of unusual phenomena, 
which have both fundamental and application interests. 
Unconventional superconductors, such as Sr2Ru04 [||, 
ZrZn 2 0] and UGe 2 @, which are characterized by the 
state with broken time-reversal symmetry (TRSB), pos- 
sess the magnetic structure related to non-zero orbital 
magnetic moment due to spin-triplet p-wave state of 
multi-component SC order parameter. The properties 
of such superconductors are not trivial. In particular, 
the macroscopic magnetization is not well-defined qual- 
ity, and the equations of motion cannot be expressed in 
terms of local magnetization [7| . 

Starting from the pioneering work [To| . [Tl| on inves- 
tigation of non-uniform states of superconducting order 
parameter (see, also review (12j). the planar defects such 
as domain walls (DW) and surfaces have been attracted 
recently a lot of interest [l3j], [HI, [H| due to new 
experiment possibilities, e. g. scanning superconduct- 
ing quantum interference device (SQUID) and Hall probe 
technique, on investigation of magnetic field distribution 
resulting from spontaneously generated superconducting 
currents which can serve as a prove of TRSB origin of su- 
perconducting order parameter in unconventional chiral 
superconductors. 

Although, the magnetostatic fields are screened by su- 
perconducting current, metastable domain walls (DWs), 
as topologically stable planar defects, may exist even in 
the Meissner state 0. Thus, the domain structure of 
SCFM cannot be ignored. The anomalies in the local 
magnetization loop for Sr2Ru04 near B = are con- 
sidered as a strong indication of the presence of chiral 



domains [15| . The strong coordinate dependence of mag- 
netization of such a 2D magnetic defect creates the intrin- 
sic magnetic field which interacts with superconducting 
(SC) current. For unconventional superconductors the 
discontinuity of magnetic induction at the DW can be 
interpreted as an effective magnetization, and the con- 
tribution of DW current to the energy density can be 
transformed to an effective Zeeman term providing the 
possibility of excitation of orbital magnetization waves 
by incident electromagnetic field 

In present publication we solve the Maxwell equations 
for distribution of magnetic field inside the sample and 
above it due to the presence of planar defect, such as 
a DW, show rapid decrease of the magnetic field in the 
vicinity of sample surface, discuss the long-range origin of 
this field resulting in existence of the tails at the distance 
~ A (London penetration depth) and compare our results 
with recent experimental data on measurements of the 
magnetic field in Sr 2 Ru04 by SQUID [l3[ and Hall probe 
technique [lj| . It is shown that for infinitely thin DW the 
stray fields above the sample depend only the magnetic 
field jump at the DW, but not on the details of DW 
structure. 



II. BASIC EQUATIONS 

Let us consider the 180° DW of Bloch type in a semi- 
infinite sample of superconducting ferromagnet occupy- 
ing V < 0. We assume, that the surface of a crystal is 
parallel to the x — z plane at y = 0, and the domain wall, 
being parallel to the y — z plane, separates two domains 
with magnetization M(x) along the +y or — y direction 
provided that the quality parameter a = Hk /AttMq > 1 , 
e.g. the field of magnetic crystallographic anisotropy of 
the easy-axis ( y-axis ) type (Hk) is high enough in or- 
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FIG. 1: Shown are the y-component (a) and z-component (6) 
of magnetic induction for Bloch DW at constant value of the 
domain wall width A and variable London penetration depth 
A. The indexes 1, 2 and 3 correspond to S — A/A = 0.1, 0.3 
and 0.6 respectively. The dashed line shows the j/-component 
(a) and z-component (6) of magnetization which are the lim- 
iting values for corresponding components of magnetic induc- 
tion at A — > co, e.g. in neglecting of screening by supercon- 
ducting current. 



der to stabilize the system against the flip of the domain 
magnetization into the x-z plane. We restrict ourselves 
to the simplest case when the London penetration length 
A exceeds the DW thickness A. This means that at the 
spatial scales of the order of A, the domain wall struc- 
ture is governed by the exchange and magnetic crystallo- 
graphic energies and is not affected by the kinetic energy 
associated with superconducting currents. Thus, we as- 
sumed that the next hierarchy of characteristic lengths 
for SCFM is fulfilled: A < A < D, where D is the size of 
the domain. We consider the steady-state properties of 
DW deep in the bulk. 

The free energy of the superconducting ferromagnet 
can be written as follows M: 



F = Jd 3 x l^na [Ml + A 2 (d i M i )(d i M i )] + ^ 



l £ V ^- A 



8^ 



B M 



(1) 



where <f> is the phase of the superconducting order pa- 
rameter, $o = hc/2e is the magnetic flux quantum, M 
is the magnetization, and Mj_ is the magnetization com- 



ponent perpendicular to the easy axis z, which is the 
normal to the sample surface. The vector potential A 
determines the magnetic induction B = curl A. The 
first and second terms in Eq. ([T]) describe the magnetic 
crystallographic and exchange energies correspondingly, 
the third term is kinetic energy of superconducting cur- 
rent, and the last two terms relate to magnetostatic en- 
ergy. The length A characterizes the stiffness of the 
spin system and is also the DW width. The magne- 
tization can be expressed in polar coordinates, M = 
Mo (sin (j> sin 9, cos 8, cos <j) sin0) with M = g^ss/a 3 the 
saturation magnetization and a the lattice constant. 



Ground state of 180° Bloch domain wall 



The ground state of the SCFM is defined by: 



5£ 

Sqi 



0, 



(2) 



where qi = 9, (j>, A. Minimization of the free energy den- 
sity Eq. |T]) with respect to azimuthal, (j>, polar, 9, an- 
gles and vector potential A leads to the system of non- 
linear differential equations describing the ground state 
of SCFM with DW of Bloch type. In case of 5 = A/A < 1 
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FIG. 2: Shown are the y-component (a) and z-component 
(6) of magnetic field for Bloch DW at constant value of the 
domain wall width A and variable London penetration depth 
A. The indexes 1, 2 and 3 correspond to 8 — A/A = 0.1, 0.3 
and 0.6 consequently. 
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FIG. 3: Shown are the y-components of magnetic induction 
(a) and magnetic field (b) for Bloch DW at constant value of 
the London penetration depth A and variable DW width A. 
The indexes 1, 2 and 3 correspond to S = A/A = 0.1, 0.3 and 
0.6 consequently. The dash lines correspond the limit values 
for magnetic induction (o) and magnetic field (6) in case of 
infinitely thin DW 0. 



the solution of these equations are as follows: 



=cos6» e (a;/A)coseQ 



6»o = 2tan~ 1 e a;/A , 

,(x/X) cos #o 

1 



00 = 0, 
=sine e 



(3) 



where we used the reduced values for the components 



of magnetic induction, b 

?(0) 



(0) 



B^/AttMq and b 



(0) 



£?z ' /AirMo. The condition 5 < 1 allows us to neglect the 
influence of Meissner current on the structure of the DW, 
which is defined entirely by the magnetic anisotropy and 
exchange energy. The screening action of the current re- 
sults in decreasing of magnetic induction at the distance 
of ~ A from the center of the domain wall x = 0. The 
graphical representation of the solutions §3§ are shown in 
Fig. [1] for different values of London penetration depth A 
at constant DW width A. The distribution of magnetic 
field H = B — 47rM, created by Meissner current, are 
shown in Fig. [2] It follows from Fig. Q] and Fig. [2] that 
the screening effect of superconducting current degreases 
with increasing of London penetration depth A. 

At constant A in limit case when the DW width A — ► 
0, the jump of the tangential component of magnetiza- 
tion My at the plane of the geometric domain boundary 
y — z defines the current sheet responsible for the jump 



of tangential component of magnetic induction B y [8j. 
The distributions of tangential components of magnetic 
induction B y and magnetic field H y of Meissner current 
at different values of the DW width A and constant Lon- 
don penetration depth A are shown in Fig. [3] The results 
of our calculation based on Eq. show the decrease of 
maximum of magnetic induction splitting with increase 
of DW width A (see, Fig. which confirms the effect 
of smoothing due to finite DW width [2(| ■ 

It follows from Maxwell equation curl H = j that 
there exists a two-component screening current }(x) — 
(0, j y (x), j z (x)), which is the result of 2-dimcnsional 
structure of the DW Eq. Q. 



III. RESULTS AND DISCUSSION 

To calculate the distribution of magnetic field near the 
surface of SCFM we neglect the domain wall A, assuming 
that A << A. This assumption does not affect the results 
qualitatively, but essentially simplifies the problem. In 
the end we shall discuss the effects of finite DW width 
A + 0. 



A. Infinitely thin DW (A -> 0) 

The difference between a normal FM and a SCFM is 
important at distances larger than A: while in normal 
FMs the magnetic field H = B 47rM vanishes and the 
magnetic induction B = 47rM is constant inside domains, 
in SCFMs the magnetic induction B is confined in the 
Meisssner layers of width A [8| : 



1 



B, = ±-8B v e ±x/x , 



(4) 



where the upper and the lower signs correspond to x < 
and x > respectively, and SB y is the magnetic induction 
splitting. Thus the Meissner currents j z = — [cj ArK)d x B v 
screen out the main bulk of domains from the magnetic 
induction. This screened magnetic induction Eq. ((4]) in- 
fluences at the distribution of magnetic field above the 
sample surface and, in principle, can be detected by 
Hall probes. To find the distribution of magnetic field 
above the surface of the sample due to magnetic induc- 
tion Eq. (|3]) the standard procedure of solving Maxwell 
equation in magnetic media (y < 0) and vacuum (y > 0) 
with appropriate boundary conditions at the surface lo- 
cated at y = 0, e. g. the continuity of normal compo- 
nent of magnetic induction and tan gen tial component of 
magnetic field, is used (see, e. g. [la]). The results of 
calculation can be represented in explicit form: 



h x = 



+°° „ fcexp 
ak — ^ — 

i) k(k + k) 



ky) 



cos (kx) , 



2 

7T Jo 



+ 00 



dk 



k exp (—ky) 
k(k + k) 



sin (kx) 



(5) 
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FIG. 4: (color online) Shown are the j/-component of mag- 
netic field in vacuum above the surface of SCFM (blue line) 
and in the sample (red line) for one Bloch DW lying in y — z 
plane. The presence of the identical tails in the distribution 
at large distances are due to long-range origin of dipole-dipole 
interaction in accordance to Eq. @. 



for the components of the reduced magnetic field hi 
2Hi/6Bi in vacuum, and 



+oc fcexp (ky) 

ak — ^7— cos ykx) , 

k(k + k) 



. 2 

Oy = — 
U 7T 



+ °° „ k[k + k(l-exp(ky))] . „ N 

dk 1 ~ J-± y ' n sin (kx) (6) 

k 2 {k + k) 



for the components of the reduced magnetic induction 
hi = 2Bi/5Bi in the sample. In Eqs. (O,© we used the 
notation k = (k 2 + A -2 ) 1 / 2 . The expression for h y (last 
equation in j5j) in slightly different form was calculated 
in 0. 

The components of magnetic field in vacuum Eqs. ([5]) 
can be calculated by introducing the complex potential 
if) = tp(w) which is analytical function of complex vari- 
able w = y — ix. Thus, the complex magnetic field 



h(w) = h y (x, y) + ih x (x, y) 



(7) 



is derived from ip(w) by h(w) = —d w ip(w). The complex 
potential tf>(w) can be expressed through special func- 
tions 



2 1 

tp(w) = -i id w [H. Q (w) 

7T W 



N (w)}, (8) 



where Hq(w) and N (w) are zero-order Struve and Neu- 
mann functions correspondingly (see, e. g. 21]). For 



\ w \> A, the difference Ho(iu) — Nq(w) w 2/ttw, there- 
fore, the main contribution to the potential tl>(w) is due 
to the first term in Eq. (jHJ which allows to calculate the 
asymptotic distribution of the component of magnetic 
field above the sample (y > A) 



-A 



y- 



xy 



7r (x 1 + y 2 ) 2 



7r (x 2 + y 2 ) 2 



(9) 



In particular, it follows from Eq. ©, that at | w \» 1 
the value of magnetic field decreases as inverse power 
law, e. g. | h \^\ w | ~ 2 , due to long-range magnetostatic 
interaction. 

The results of numerical calculations of normal to the 
sample surface component of magnetic field in vacuum 
above the surface of the sample h y Eq. © and in the 
sample b y Eq. (JSJ) are represented in Fig. 3] which shows 
that at the sample surface y =0, magnetic induction split- 
ting equals half of it bulk value SB y and are characterizes 
by the rapid decrease with the distance above the DW. 
Thus, at the distance of y =0.1 A above the DW the mag- 
netic field decays till the third of the bulk value 6B y . At 
y ~ A, the magnetic field splitting reaches about 0.05 of it 
bulk value and decays with x as inverse power low Eq. ^ 
due to the long-range origin of dipole-dipole interaction. 

In a limit of infinitely thin DW, A — > 0, the stray 
fields in vacuum above the superconducting ferromagnct 
depend only on the jump of magnetic induction at DW 
plane, but not on the details of DW structure. In this 
case, the tangential to the sample surface component of 
magnetic field is non-analytic function of y and is char- 
acterized by the discontinuous jump at the DW position, 
y = 0. In next subsection we take into account the effects 
of small (A << A), but finite DW width A ^ 0. 



B. Finite DW width (A / 0) 

To consider the effects related to the finite DW width, 
A / 0, we assume the linear distribution of magnetic 
induction in the DW deep in the bulk, e.g. B y w 
(l/2)8B y (x/ A) at | x |< A, which is not affected by 
superconducting currents, inasmuch as A << A. This 
assumption significantly simplifies the calculation and al- 
lows to express the results in analytic form. It can be 
shown that for finite DW width the kernels in Eqs. ([5]) 
and ([6]) are modified by the multiplier sin(fcA)/fcA. It 
results in slight suppression of y-component of magnetic 
field and restores the analyticity of x-component of mag- 
netic field. The discontinuous jump of tangential compo- 
nent of magnetic field at the position of DW is replaced 
by the value h™ ax = (1/tt) ln(jS), where 7 « 1.78107 is 
Euler-Mascheroni constant. 

Recent experiments on imaging of magnetic field dis- 
tribution above the surface of the sample of S^RuC^ 
by scanning SQUID and Hall probe microscopy has re- 
vealed no evidence for existence of DWs in this uncon- 
ventional superconductor (l9j . This negative result can 
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be understood in the framework of the developed the- 
ory. If the DWs exist, the maximum value of magnetic 
induction due to permanent current flowing in the plane 
of the DW does not exceed the lower critical field, e. g. 
SBy/2 ~ H c i «30G above which the domain structure is 
unstable due to formation of Abrikosov vortices. Thus, 
at the distance of y = A «190 nm, the normal component 
of magnetic field due to the presence of DW is of order 
of 3G which definitely improbable by the used methods. 
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